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Introduction 

The objectives of this project are to gain insight into diffusion flames by modeling various 
configurations related to ongoing experimental investigations in the microgravity combustion 
science program. The emphasis of the work is to understand the structure and dynamics 
of diffusion flames. Improving our fundamental understanding of diffusion flames is most 
relevant to issues related to fire safety and fire prevention because most fires consist of diffusion 
flames. 

Work on this project started in May 1994. 


Npar-limit Oscillations 

From recent experiments of candle flames conducted in a reduced gravity environment 
[1], [2], the following observations were made: 

• the flame shape is nearly a hemisphere as opposed to the elongated tear drop shape 
observed under normal gravity conditions 


• the flame standoff distance is quite large; at the base of the flame that distance is about 
5 mm compared to 1 mm in normal gravity conditions 


• the flame is steady during most of the process, however, as the oxygen in the chamber 
is cons um ed, spontaneous oscillations develop prior to extinction 

• the burning rate and the flame temperature are considerably reduced when compared 
with their values in normal gravity. 


The goal of our current efforts is to shed light on these observations and in particular on the 
oscillations that occur during the last stages of the burning. 

In microgravity the candle flame bears some similarity to the diffusion flame that sur- 
rounds a burning liquid droplet. Both have a spherical shape and both depend on the heat 
conducted back to the condensed phase. Hence their structure and dynamical properties are 
likely to be similar. In our first attempt to model the above observations we have therefore 
analyzed the droplet problem. To simulate the experimental conditions we consider the burn- 
ing to occur in a reduced oxidant environment, i.e. the ambient mass fraction of the oxidant, 
Xoo, is assumed small. Under these conditions the flame standoff distance is relatively large 
and the burning rate is reduced. The combustion field is now divided into a near-field, where 
fuel vapor in the absence of oxidant is transported outward by convection and diffusion, and 
a broad chemically reactive far-field where convection is effectively negligible. For steady 
burning the mass burning rate rh is expressed in the form 
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where a is the droplet radius, A and c v are the thermal conductivity and specific heat (at 
constant pressure) of the mixture, B = c p (T 00 - T s )/L is the transfer number for pure 
vaporization with T 3 ,T <X> the surface and ambient temperatures respectively, L is the latent 
heat of vaporization and Q / v\ Wjf is the heat of reaction per unit mass of oxidant consumed 
(the subscript X identifies the oxidant). The parameter n, obtained by numerically solving 
the reactive-diffusive far-field equations, depends on the physico-chemical parameters and in 
particular on the Damkohler number 6 representing the ratio of diffusion time to chemical 
reaction time. We note that in the absence of heat loss fj, — *• 1 as 6 —* oo; eq. (1) then 
reduces to the classical expression for the burning rate (e.g. [3]) when expanded for small 
Xoo. The dependence of the mass burning rate n on the Damkohler number 6 is shown in Fig. 
1. The response curve exhibits different characteristics depending on the ambient oxidant 
concentration X ^ and on the heat loss parameter h. For moderate values of X^ the familiar 
S-shaped response curve is obtained. However, by reducing the oxidant concentration the 
response curve becomes monotonic connecting the weakly burning state where n ~ <5 1/2 (for 
small 5), to the equilibrium limit where fi tends to a constant < 1 as 6 — *• oo. The asymptotic 
value reached for large 6 depends of course on the heat loss parameter h. The steady problem 
in the absence of heat loss and with a unity Lewis number has been previously discussed in 
[4], as part of a complete discussion of a droplet combustion analysis. 

The change in shape of the response curve has significant implications regarding the 
dynamics of the flame. For example, the concept of abrupt extinction is associated with 
the S-shaped response curve. If the Damkohler number 6 is slowly decreased, say in an 
experiment, the upper branch corresponding to an intense burning is followed down to the 
turning point Se where the response jumps to the lower branch corresponding to a weak 
burning or an extinguished state. The point 6 = 6e (which, in general, does not have to 
be the turning point) is associated with extinction. In contrast, a monotonic response curve 
indicates that the transition from an intense to a weak burning takes place in a gradual way. 
Behind this description is the assumption that the middle branch of the S-curve is unstable 
and therefore cannot be realized in practice. The monotonic curve on the other hand is often 
assumed to be stable for all 6. This hypothesis appears to be true when there is no heat loss. 

In the presence of heat loss, however, it is found that part of the monotonic curve may 
become unstable (see Fig. 2). When perturbations proportional to exp(iu;t) are sought, one 
finds that growing modes are possible for 6 below some critical number. Furthermore, at 
the onset of instability Ke(u) ^ 0 implying that a growing oscillatory mode develops. The 
frequency of oscillations is given by 


7 Ze{u>) = a 
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where (3 = E/RToo is the activation energy parameter, assumed large, and a is an 0(1) 
number. By slowly decreasing 6 below its critical value, only a small change in the frequency 
(i.e. in o) is detected but the amplitude grows continuously. It is therefore anticipated that as 
the amplitude becomes sufficiently large the flame extinguishes. Hence near limit oscillations, 
similar to the one observed in the microgravity candle flame experiment [2], may result in 
the presence of heat loss. 

Future plans include examining other source of perturbations which may also yield flame 
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oscillations. For example, near-limit oscillations have been observed in candle flames under 
normal gravity conditions at relatively low pressures [5]. It has been suggested that the 
presence of convective flows may be responsible for these oscillations. Although this does not 
appear to be the mechanism responsible for the oscillations in the microgravity experiments 
[2], it would be of interest to incorporate weak buoyant flows as a source of perturbations 
in the analysis in order to better understand the difference between these two environments. 
Attempt will also be made to model the candle flame itself, which more closely resembles a 
spherical cap rather than a complete sphere. 

Flame Front Instabilities 

Studies of intrinsic instabilities in flames have been predominantly concerned with pre- 
mixed systems. A manifestation of an intrinsic instability is the spontaneous development of 
cellular structures, which is commonly observed in premixed flames. The competing effects 
of thermal and mass diffusivities play a central role in the development of cellular flames. 
One might therefore expect that they play a similar role in diffusion flames as well. However, 
very little has been done on the subject. The lack of experimental work in this area may be 
attributed to the fact that it is difficult to generate in the laboratory a one-dimensional diffu- 
sion flame free of the stabilizing influence of stretch. The most comprehensive experimental 
study [6] has been carried out only recently using a slot- jet flame. Using a spherical burner 
in microgravity a one dimensional spherical flame of sufficiently large size can be produced. 
This configuration may well be suited for studying the development of flame front instabilities 
in diffusion flames, both theoretically and experimentally. The goal of our current efforts is 
to pursue a theoretical study of thermal-diffusive instabilities in diffusion flames. 

An important parameter in this study is the Damkohler number 8. In the limit 8 -*• oo, 
known as the Burke Schumann limit, the chemical reaction is concentrated in a thin sheet 
separating a region where there is fuel but no oxidant from a region where there is oxidant 
but no fuel. The requirement that fuel and oxidant do not co-exist and the requirement that 
fuel and oxidant must flow into the flame sheet from opposite sides in stochiometric propor- 
tions, are sufficient conditions for complete determination of the combustion field including 
the instantaneous shape and location of the flame sheet. Stability analyses of simple flame 
configurations show that under such conditions the flame sheet is absolutely stable. 

The situation is different when there is a leakage of fuel and/or oxidant. The mathe- 
matical formulation in this case is based on the assumption that 8 is large but such that 
6 = 8/3~ 3 exp(-/3/Tf), where /? > 1 is the activation energy parameter and T/ is the flame 
temperature for complete combustion, is 0(1). The reduced Damkohler number 8 plays an 
important role in this description. The inner structure for a steady one-dimensional flame 
has been discussed by Linan [7]; it depends on two parameters 8 and 7. (Of practical interest 
is I7I < 1; depending on its sign 7 stands for the fractional excess of heat that is conducted to 
the fuel or oxidant side of the reaction zone). Numerical integration of the structure equation 
reveals the dependence of the combustion field on these parameters. The Burke Shumann 
limit is recovered when 8 -*• 00. By reducing 8 a small leakage of fuel and/or oxidant develops 
and the flame temperature decreases. This behavior persists with the reactants leakage and 
the flame temperature reduction increasing continuously, until extinction occurs at 8 = 8 ex . 

For a steady one-dimensional flame the parameter 8 is a constant determined by the 
underlying flow field. For example, in a counterflow 6 is proportional to the imposed strain 
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rate. For a general three-dimensional corrugated flame 6 is only determined locally, for 
example in terms of the local orientation of the reaction sheet and the local fuel and oxidant 
concentrations, or in terms of the local gradients of the mixture fraction [3]. We have thus 
extended Linan’s analysis to allow for temporal and spatial variations. To complete the 
mathematical formulation expressions for the reactant leakages and their dependence on 6 
must be identified. Fig. 3 shows the results of the numerical integration of the structure 
equation for 7 > 0. For 7 < 0 the role of the fuel and oxidant are interchanged. Expressions 
for the leakages that best interpolate the numerical data of Fig. 3 take the form 

Y. ~ (PS)' 1 exp - ^ei)°] (2) 

X , ~ /3 -1 exp - £ e x) fc ] ( 3 ) 

where CF,Cx,a,b depend only on 7. Clearly Y. and X. vanish as 6 -> 00. When 7 = 0, 
equal amounts of fuel and oxidant leak through the reaction zone. As 7 — ► 1, however, the 
fuel leakage at extinction becomes excessively large, i.e. Y, = 0(1), and the oxidant leakage 
becomes negligibly small, i.e. X * = 0. The reverse is true when 7 — *• — 1. This general formu- 
lation is suitable for investigating the stability of simple flame configurations. Preliminary 
results indicate that depending on the fuel and oxidant Lewis numbers instabilities leading 
to cellular structures may arise. This appears to be consistent with the experimental results 
of [6]. A systematic study is being carried out to identify the stability boundaries in terms 
of the relevant parameters. 
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Fig.l - The dependence of burning rate on the Damkohler number 
for selected values of 7 ~ Xx,; calculated for h = 0. 



Fig.2 - Response curves for selected values of 7 ~ Xoo with h — 1. 

The dark segments correspond to the unstable states where the oscillations occur. 
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